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ABSTRACT 

We have detected the tip of the Red Giant Branch (TRGB) in the solar neighborhood using near 
infrared photometry from the 2MASS and DIRBE catalogs, and revised Hipparcos parallaxes. We 
confirm that the revised Hipparcos parallaxes are superior to the original ones, and that this improve- 
ment is necessary to detect the TRGB. We find a tip absolute magnitude of Mk = —6.85 ± 0.03, in 
agreement with that expected from previous tip measurements of the Large Magellanic Cloud, Small 
Magellanic Cloud, and Bulge. This represents the first geometric calibration of the TRGB and extends 
previous calibrations, based on metal-poor globular clusters, to solar metallicities. We attempted to 
use the TRGB to confirm the presence of the Lutz-Kelker bias, with inconclusive results. Attempts to 
detect the tip in the /-band also produced inconsistent results, due to a lack of precise, homogeneous 
photometry for these bright stars. 

Subject headings: stars: AGB and post-AGB - stars: late-type - stars: variables: other - solar 
neighborhood - distance scale 



I. INTRODUCTION 

The tip of the Red Giant Branch represents the max- 
imum absolute luminosity achieved by first-ascent red 
giants, and marks the onset of helium fusion in the 
degenerate cores of these low-mass stars. Theoretical 
predictions, confirmed by observational evidence, indi- 
cate that the TRGB is an excellent distance indicator, 
since the absolute bolometric magnitude of the tip varies 
by only ~ 0.1 ma g for a wi de range of metallicities 
and ages (llben fc R enzini 1 983tlDa Costa fc Armandrofil 
fl990l: ISaJaris fc Cassisitd997ir The /-band magnitude of 
the TRGB has become well-established as a distance in- 
dicator for nearby galaxies with well-res olved Population 
II stars in their halos (for example, se e iLee et all 119931 : 
ISakai et al.lll996l : IMould fc Sakaiil2008h . 

The TRGB complements other distance indicators, 
such as Cepheid s and RR Lyra e stars, being compara- 
ble in precision ()Lee et al.l ll993'). Indeed, empirical evi- 
dence and computer simulations show that space-based 
observations provide distances out to ~12Mpc with a 
precision of 10%, li mited primarily by integration time 
IjMadore fc Freedm an 1995). However, the TRGB is a 
tertiary distance indicator, calibrated using RR Lyrae 
variables in Galactic globular clusters, which are them- 
selves secondary distance indicators. Moreover, there is a 
discrepancy of ^O.lmag between the T RGB and Cepheid 
distance scales (iTammann et al.|l2008[ ). and the metallic- 
ity dependence of Cepheid P-L relation has itself been 
calibrated using the TRGB meth od, leading to a circu- 
lar dependency (|Rizzi et all 120071 ). Thus, it would be of 
considerable importance to have a direct calibration of 
the tip magnitude. 

Although the TRGB is more easily observed in near in- 
frared (NIR) passbands, where interstellar reddening is 
reduced, both the color and luminosity of the NIR TRGB 
are more sensitiv e to metallicity than in the /-band 

izzi et ahl |2007| ). Nevertheless, recent studies of the 
Magellanic Clouds and Galactic Bulge have successfully 
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identified the TRGB discontinuity in the /T-band, using 
2MA SS photometrv (ICioni et al.l l2000l: iKiss fc Beddh^ 
12004 ISchultheis et al.ll2004f ). Moreover, the discontinu- 
ity is clearly visible in the /tT-b and P-L plot of M g iants in 
the local solar neighborhood (T abur et al.l 120091 ) . which 
prompted us to investigate further. In this Letter we 
present our findings, including the first geometric cali- 
bration of the TRGB absolute magnitude. 

2. DATA 

We selected all stars with relative parallax uncer- 
tainties less than 25 % in the revised Hipparcos catalog 
(jvan LeeuwenI l200"7i )^ . We obtained JHK mag nitudes 
and u ncertainties from the 2MASS catalog (Cutr ret al.l 
[2003h . using a search radius of 5 arcsec. Being luminous 
and nearby, some of the stars were saturated despite be- 
ing observed with the shortest integration time (51ms), 
and were measured using the wings of their radial pro- 
files, resulting in large photometric uncertainties (0.2- 
0.4 mag). Additional scatter was contributed by the use 
of single-epoch mea surements, since all M giants are in- 
trinsically variable (jEver fc Mowlavill2008l ). with typical 
/C-band peak-to-peak amplitudes between 0.1 and 0.25 
mag (Smith et al. 2004). 

Seeking a higher-precision source of NIR photometry, 
we extracted fiux measurements from the DIRBE cat- 
alog, which contains photomet ry for nearly 120 00 ob- 
jects sampled over 10 months (jSmith et al.fl2004D . Fol- 
lowing Whitelo ck et al.l ()2008f ). we adopted a value of 
630 Jy for niK — 0, and used 1570 Jy for mj = , 
both of which are consistent with lBessell fc BrettI (|1988f ). 
Although DIRBE was less prone to saturation than 
2MASS, its large beam-width led to confusion between 
nearby so urces. Using the same methodology as de- 
scribed bv lTabur et al.l (|2009f) . we only selected uncon- 
fused sources (confusion flags 1 and 2 not set). This per- 
mitted an accurate determination of mean K magnitudes 

^ We used the version published on VizieR on 2008 September 
15, which corrected an error that affected earlier versions of the 
revised catalog. Stars with negative parallaxes were ignored. 
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that are demonstrably superior to 2MASS for about half 
the red giants. 

Interstellar extinction is expected to be relatively small 
for all stars in the sample, particularly in the K-hand. 
We calculated the visua l extinctions us i ng the interstel- 
lar extinction model of iDrimmel et all (|2003f ). Extinc- 
tions for other bands we re scaled using the factors from 
iRieke fc Lebofskvl ()1985[ ). The greatest extinction was 
Ak ~ 0.15 mag, with Ak < 0.05 mag for 98% of the 
sample. 

We calculated absolute magnitudes using the relation 
M = m + 5 + 5 logTT, where m is the extinction-corrected 
apparent magnitude, and tt is the geometric parallax in 
arcseconds. Photometric uncertainties were obtained di- 
rectly from the 2MASS catalog, or derived from flux 
densities from DIRBE. Parallax uncertainties were esti- 
mated with a first-order approximation of 5 loge ct^/tt w 
2.17 cttt/tt, with the total uncertainty in M calculated as 
the quadrature sum of photometric and parallax uncer- 
tainties^. 

The color-magnitude diagrams (CMDs) are shown in 
the top row of Figure [U and the RGB is clearly vis- 
ible. We indicate its expected position with a thin 
black line, by assuming [Fe / H] = and adopting the 
calibration of IValenti et al.l (|2004f ). which is based on 
Galactic globular clusters within the metallicity range 
-2.12 < [Fe/H] < -0.49. The position agrees weh for 
our 2MASS sample, and less so for DIRBE, although the 
difference is irrelevant for measuring the TRGB. A few 
stars around J ~ K ^ 1.5 and Mk ~ are probably 
reddened by dust. The DIRBE photometry is sparser, 
particularly at the faint end, although the evolved stars 
near the tip appear well represented. 

To be sure of identifying all stars near the tip, which 
contains relatively few stars because it is a short-lived 
late evolutionary stage, we ignored cataloged spectral 
types. Instead, we used the color-magnitude diagram 
to select all stars in the range 0.5 < J — K < 2.0 and 
Mk < —3, in order to include early M giants, and to 
exclude very red, dust-enshrouded stars. 



3. EDGE DETECTION 

We have used two model-independent methods for 
edge detection (ED) to l ocate the TRGB d iscontinu- 
ity. Firstly, we followed iMager et all (|2008l hereafter 
MMF08) by using an ED filter that uses a combination 
of small bins for resolution and larger bins for smooth- 
ing, with a logarithmic edge detector to account for the 
expect ed power-law distrib ution at the bright end of the 
RGB (|M6ndez et al.ll2002[) . Additionally, we iteratively 
calculated the peak filter-response 85 times, using a range 
of starting offsets and bin sizes to smooth-out sensitiv- 
ity to these parameters, and subsequently used the mean 
(after 3-sigma rejection) as the TRGB location. This 
method proved fairly robust in the presence of noise. 
We limited the search for the discontinuity to the range 
— 8 < Mk < —6, to eliminate false positives. 

Our second method accounts for photometric and par- 
allax uncertainties by using a sum of Gaussians to create 

^ The data used in this analysis are available electronically from 
CDS, Strasbourg. 
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Fig. 1. — Top: Color-magnitude diagrams using photometry 
from 2MASS (left) and DIRBE (right) for all stars with revised 
Hipparcos parallaxes satisfying (t^/tt < 0.125. Middle: luminosity 
functions for red giants, with bin sizes as shown. Bottom: results 
of applying two edge-detection methods to the LPs (see text). 

continuous probability function to model the LF: 
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where and ai are the magnitudes and their uncer- 
tainties, respec tively, and N is t he total number of stars 
in the sample (|Sakai et al.lll996l hereafter SMF96). We 
found that a logarithmic ED filter produced poor re- 
sults and, following SMF96, adopted a filter of the form 
E{m) = ^(m -|- Am) — <f>(m — Am), where Eirn) is the 
filter response at magnitude to, and Am is the bin size. 
Previous studies have defined the bin size as a function of 
photometric uncertainty. However, unlike the RGBs in 
nearby galaxies, our sample is drawn from the local solar 
neighborhood and is less affected by completeness errors 
or increased photometric errors toward fainter magni- 
tudes. Indeed, the distribution of mean errors over the 
range —8 < Mk < —6 is nearly constant, leaving us 
free to choose a bin size on the basis of its smoothing 
performance alone. We adopted a value of 0.2 mag. 

We estimated unc ertainties using the boo tstrapping 
method described bv lBabu &: Feieelsonl (Il996|). wh ich is 
widely used (see, for example. iMendez et al.ll200l) . We 
selected 1000 random samples, each consisting of 80% of 
the stars from the original sample, and determined the 
location of the TRGB for each, adopting the rms scatter 
as our uncertainty. 
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TABLE 1 

X-BAND TRGB MAGNITUDES FROM 2MASS AND DIRBE. 



Max 


N2M 




. ,MMl-',i!M 




, .MMV.D 


CTtt /tT 




(mag) 


(mag) 


(mag) 


(mag) 


0.100 


123 


-6.83 ± 0.03 


-6.79 ± 0.02 


-6.82 ± 0.03 


-6.84 ± 0.04 


0.125 


166 


-6.83 ± 0.03 


-6.82 ± 0.03 


-6.82 ± 0.03 


-6.84 ± 0.04 


0.150 


212 


-6.85 ± 0.05 


-6.87 ± 0.03 


-6.82 ± 0.03 


-6.90 ± 0.09 


0.175 


253 


-6.85 ± 0.03 


-6.85 ± 0.03 


-6.82 ± 0.03 


-6.85 ± 0.06 


0.200 


295 


-6.86 ± 0.03 


-6.83 ± 0.03 


-6.82 ± 0.03 


-6.85 ± 0.06 


0.225 


318 


-6.86 ± 0.03 


-6.82 ± 0.03 


-6.82 ± 0.03 


-6.86 ± 0.06 


0.250 


353 


-6.87 ± 0.03 


-6.85 ± 0.05 


-6.82 ± 0.02 


-6.94 ± 0.10 



TABLE 2 

if-BAND TRGB MAGNITUDES CORRECTED FOR LuTZ-KeLKER BIAS. 

Max Al^ ' ' ' A/^ 

o-j/tt (mag) (mag) (mag) (mag) 

0.100 -6.86 ± 0.03 -6.81 ± 0.02 -6.84 ± 0.03 -6.86 ± 0.02 

0.125 -6.86 ± 0.03 -6.83 ± 0.03 -6.84 ± 0.03 -6.86 ± 0.03 

0.150 -6.87 ± 0.03 -6.93 ± 0.05 -6.83 ± 0.04 -6.86 ± 0.05 

0.175 -6.86 ± 0.03 -6.84 ± 0.04 -6.83 ± 0.03 -6.82 ± 0.03 

0.200 -6.87 ± 0.03 -6.85 ± 0.04 -6.83 ± 0.03 -6.83 ± 0.02 

0.225 -6.87 ± 0.02 -6.84 ± 0.03 -6.83 ± 0.03 -6.83 ± 0.02 

0.250 -6.89 ± 0.02 -6.83 ± 0.03 -6.83 ± 0.03 -6.81 ± 0.02 

(|2000f ) measured the location of the TRGB in 10 galac- 
tic globular clusters, fitting the relation M™*^^ = 
- (0 . 59 ± . 1 1 ) [Fe/H] ^^^^ - (6.97±0.15), which predicts 
mag for an increase in metallicity by 
a factor of 2. Using the global me tallicity scale, which 
incorporates a-element abundances. iFerraro et al.l ()2006D 
found M^i'GB ^ _g 92 _ o.62[M/H] for metalhcities to 
[Fc/H] < +0.4 dex. This leads to a tip magnitude in the 
range —6.8 to —6.9 for solar-like metallicity, in excellent 
a greement with ou r resul t . 

ISchultheis et al.l (|2004f ) examined the K-hand lumi- 
nosity function of MACHO semiregular variables in the 
Galactic bulge, finding a steep drop-off in number near 
Mk ^ —6.8, which they interpreted as the location of 
the TRGB. 

4.2. Lutz-Kelker bias 

We expect a sample limited by relative parallax 
uncertainties to be affected by a systematic bias 
(jLutz_& Kelkeij [19731), resulting in an underestimate of 
stellar luminosity. Using the TRGB as a marker, we 
searched for evidence of a systematic shift in its position 
in samples selected with progressively larger relative par- 
allax uncertainties. 

First, we selected non-overlapping samples in blocks 
of 10% relative parallax uncertainty, but this yielded too 
few stars for reliable TRGB detection. Increasing to 15% 
increments did not help, so finally we simply examined 
our previous samples for a systematic shift, and com- 
pared them to a reprocessed set, where individual bias 
corrections had been applie d to each star using t he re- 
lation LK = -8.09(cr^/7r)2 (jWhitelock et al.l|2008D . Ta- 
ble O lists the results using the same nomenclature as 
Table [U and shows that the mean corrected tip posi- 
tions are indistinguishable from the uncorrected values, 
to within the errors. We note that the corrected tip mag- 
nitudes for 2MASS are 0.02 mag brighter on average, but 
the DIRBE MMF08 resuh is slightly fainter. 

The median a-jt/iT for each sample is about half the 
maximum value, implying an expected bias-induced shift 
of '--^O.l mag between the first and last samples (10% and 
25%). However, the uncorrected results (Table [1]) show 
far smaller differences, and in the wrong sense (brighter 
tip magnitudes when an under-estimate is expected). 
Thus, we see no evidence of a systematic shift toward 
fainter tip magnitudes as samples include stars with pro- 
gressively larger relative parallax uncertainties, and con- 
clude that there are too few stars with sufficiently precise 
parallaxes and/or apparent magnitudes for a convincing 
detection of the Lutz-Kelker bias. 

4.3. Comparison of original and revised parallaxes 



4. RESULTS AND DISCUSSION 
4.1. The K-band TRGB 

Since samples containing stars with relatively large 
(Ttt/tt tended to blur the RGB/AGB boundary, we cal- 
culated the TRGB location for several samples selected 
with upper bounds in the range 0.10 < aTr/n < 0.25. Lu- 
minosity functions (LFs) and bootstrapping distributions 
for the stars with (Ttt/vt < 0.125 are shown in Figure [TJ 
LFs are shown for bin sizes of 0.2-0.5 mag, together with 
the weighted continuous function (Eq. [1]). The two LFs 
exhibit similar morphologies and, despite the DIRBE LF 
containing half as many stars as 2MASS, both show the 
AGB stars above the tip, at —9 < Mk < —7. 

The bottom panels show the TRGB position deter- 
mined from 1000 bootstrap operations for each ED fil- 
ter, binned in 0.01 mag increments. The 2MASS pho- 
tometry yields an unambiguous detection, with nearly 
identical results for the two filters. The sparser DIRBE 
photometry produced very similar results, although the 
SMF96 filter detected another peak at Mk ~ —7.1 which 
was not real. Its significance reduced as more stars were 
added (with looser parallax restrictions), simultaneously 
increasing the size of the peak at Mk ^ —6.8, imply- 
ing that the SMF96 filter is more sensitive to noise. We 
ignored the second peak in our calculations. 

Table[l]gives the full results, listing the maximum (Ttt/ 
used for each sample, the number of stars in the 1-mag 
bin below the tip (2MASS), the tip locations for both 
the 2MASS and DIRBE samples using the SMF96 and 
MMF08 filters. The results for all four combinations 
agree within the uncertainties. 

Early studies of the statistical robustness of the TRGB 
method concluded that 50-100 stars are required in the 
1-mag bin below the TRGB, although more recent work 
suggests that 4 00-500 stars are req uired for an accurate 
determination (jMadore et al.ll200"9l ) . Although our small- 
est 2MASS sample includes only 123 stars in the 1-mag 
bin below the tip, the largest sample (with 353 stars) 
yields very similar results. There is no indication of the 
0.6 mag systematic shift toward fainter tip magnitudes 
for samples with less than 300 stars, that was found 
bv lMadore et al.l (|2009D . Moreover, rerunning the ED 
calculations by bootstrapping with 66% and 33% of the 
original sample size showed little variation in the TRGB 
lo cation (but with l arger uncertainties). 

ICioni et all (j200a) measured the TRGB in the LMC 
and SMC, finding dereddened values of Mk = —6.61 
and —6.41, respectively. The difference of 0.2 mag is 
caused by the different metal abundances, which we 
adopt as Zcd = 0.016, Ztmc. = 0.008, and Zsmc = 
0.004 (jVassiUadis fc WoodI Il99l . In or der to quan- 
tify the TRGB's metallicity dependence, iFerraro et al.l 
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Fig. 2. — Same as the bottom row of Fig. [T] but with Mk 
calcu lated using the original Hipparcos parallaxes (Perryman et al. 
119971) . As in Fig. [T] the left and right panels are for 2MASS and 
DIRBE photometry, respectively. 

To test whether the revised Hipparcos parallaxes are 
superior to the original values, we use the ability to 
resolve the TRGB, with consistent results, as an in- 
dicator of precision. We reprocessed the same sam- 
ple used for Figure [H but using the original parallaxes 
(jPerrvman et alJ Il997t l. Figure [2] shows the resulting 
bootstrap distributions for 2MASS and DIRBE photom- 
etry. 

We note several features. Firstly, the peaks are no 
longer coincident, with the filters converging to differ- 
ent values. Some peaks are asymmetrical, with long 
tails, and in the case of SMF96 with DIRBE data, mul- 
tiple solutions are found. None of the tip locations 
match the value consistently returned using revised par- 
allaxes. Moreover, the peaks have significantly larger 
widths (0.05-0.12 mag, excluding the worst case with 
multiple solutions). Comparing with Figure [1] (bottom 
row), and noting the change in scale, we see that only the 
revised parallaxes produce a consistent result with small 
uncertainties. Tests using larger samples with greater 
relative parallax limits produced similar results. We con- 
clude that the revised Hipparcos parallaxes are superior 
to the originally published values, and that this improve- 
ment is necessary for a clear detection of the TRGB. 

4.4. The I-band TRGB 

An /-band calibration of the tip is particularly impor- 
tant because of its relative insensitivity to metallicity and 
age. Unfortunately, we are not aware of a homogeneous 
source of all-sky /-band photometry for these very bright 
stars. For example, DENIS photometry is saturated at 



771/ ^ 10 (|Epcht em et al.lll994f) . Instead, we transformed 
2MASS NIR magnitudes to the /-band using the rela- 
tion 777/ = 0.0560 -h 2.0812 TOj-h 0.4074 77i//- 1.4889 rriK 
(G. Bakos, pers. comm.), with uncertainties of ~0.1- 
0.2 mag. This required a tighter constraint on Ut^I-k to 
compensate for the larger photometric uncertainties in- 
troduced by the transformation, which reduced the num- 
ber of stars near the tip. 

Subsequent edge searches yielded inconsistent results. 
Samples having a-^/i: < 0.1 produced very broad boot- 
strap distributions around A// ~ —4.1 with vague, ill- 
defined peaks. Samples with greater relative parallax 
uncertainties produced slightly more consistent results, 
with Ml generally centered at ^ —3.8 ± 0.2 mag, al- 
though inspection of the corresponding LFs showed no 
obvious discontinuity. Clearly an improved source of pre- 
cise, all-sky /-band photometry is required before at ac- 
curate calibration can be made. 



5. CONCLUSION 

The TRGB is an important, widely- used tertiary dis- 
tance indicator, but lacks a direct calibration. Using re- 
vised Hipparcos parallaxes and NIR photometry from the 
2MASS and DIRBE catalogs, we have unambiguously de- 
tected the TRGB in the local solar neighborhood. Two 
different ED methods have been used to quantitatively 
measure its location as Mk = —6.85 ±0.03, which is the 
first geometric calibration of the TRGB, and the first for 
stars with solar metallicity. We demonstrate that the re- 
vised Hipparcos parallaxes yield consistent results with 
far smaller uncertainty than the original Hipparcos par- 
allaxes, and conclude that the revised values are indeed 
superior, at least for our sample of nearby, M giants. We 
were unable to detect the effect of the Lutz-Kelker bias 
using the TRGB location as an indicator. Calibration of 
the /-band TRGB is not possible due to a lack of suSi- 
ciently precise photometry. 

This research has made use of the data products from 
the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
This project has been supported by the Australian Re- 
search Council. 
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